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Abstract 

An influence of initial distribution nonequilibrium parameters on ther- 
modynamical equilibrium recovery processes in early Universe under the 
^— I assumption of elementary particles' scaling of interactions in range of su- 

■ perhigh energies is researched. 



1 Introduction 



In previous papers of the Authors the model of superheat particles' cosmolog- 
Q ' ical evolution in conditions of scaling of interactions was built up [U [2]. For 

elementary particles' cut set of interactions at that the Universal Asymptotic 
Cut Set of Scattering, introduced in papers was used: 
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00 where s is a kinematic invariant of four-piece reaction (details see in [1]), sq 

is a square of two colliding Planck masses' total energy, 

cn 

O ■ A(s) = 1 + In^ — « Const, (2) 

So 

In [5] a relaxation of cosmological plasma's superheat component on equilib- 
rium component was researching under the assumption that particles number 
in superheat component is far less than the equilibrium particles number. In 
particular, it was shown that the solution of kinetic equation describing an evolu- 
. tion of the ultrarelativistic superheat component in the equilibrium cosmological 

$H ' plasma has a form: 



AUt,V) = Af,{V)exp 



at,V) fyHt')dt' 



V J ^ 



(3) 



where Af^iV) ~ A/a (0,7-") is an initial deviation from equilibrium (detail see 
in [2]) and the dimensionless function is introduced: 

yit) = ^y ^^y'; ^o^y''(0). (4) 



^system of units G = h - 
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where T(t) is a temperature of plasma's equilibrium component, To(t) is a tem- 
perature in the same point of time in completely equilibrium Universe and the 
dimensionless momentum variable V (conformal momentum) is introduced: 



Farther M is an efhcient number of particles' equilibrium types; 
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3/ArV327r3; A{VTTq/2)' 



(5) 



(6) 



is a parameter, weakly dependent on variables t,V, 



^(25+1) + 1^(25+1) 



Nb + ^Nf; 



Nb is a number of sorts of equilibrium bosons, Np is a number of sorts of equi- 
librium fermions. Relative temperature y{t) of plasma's equilibrium component, 
presenting a parameter of the nonequilibrium distribution ^ , is determined via 
the integral equation of energy-balance: 



oo 



X exp 



at,V) f y\t')dt' 
V J ^ 



dV = 1, 



(7) 



where {2Sa -I- 1) is a statistical factor. From ([7]) in zero point of time follows the 
relation |2]: 



oo 

Y^{2Sa + l)J V^AfAV)dV = 1 - ao. 



(8) 



Let us introduce corresponding to [2] a new dimensionless time variable r : 



(9) 



where ^ = C(^o); 



E(25a + 1) / dVV'Af^iV) 



E(25„ + 1) / dVV^AfO{V) 

a 



(10) 
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is an average value of momentum variable V in point of time t = 0, Z{t)\ 

T 

Z{t) = 2 / y\T')dT'. (11) 





Then subject to ([5]) equation ([7]) is reduced to the form([5]): 

y=[l-{l-a^mZ)Y'\ (12) 
where 

oo 

= = JL 2 _ . (13) 

a 

2 Numerical Model 

In paper [5] distribution function's initial deviation from the equihbrium was 
represented in the most plain form: 

^/"(") = p3(fc2f, 2)3/2 (14) 

where xl-^) is a Heaviside function (a staircase function), x = V /Vo is a dimen- 
sionless momentum variable, A, Vo and k are certain parameters; k parameter 
is introduced to provide a convergence of distribution function's all moments in 
range of momentum's small values. A distribution function of superheat parti- 
cles' energy density at that has a form, close to the spectrum of the so-called 
white noise, when all energy's values are equiprobable up to the certain crit- 
ical value Vq, after which distribution breaks. From the qualitative point of 
view an evolution of the initial distribution is reduced to the "corrosion" 
of nonequlibrium distribution spectrum's low-energy part while energy density's 
constancy in high-energy part conserves, what did not allow to make predictions 
about the form of distribution's tail area. In given paper we select a distribution 
function's deviation in more realistic form, allowing to carry out a research of 
the relaxation process dependence on parameters of the initial distribution. The 
paper's main goal at that is a stability's research of the nonequlibrium cosmo- 
logical scenario '2] with respect to the parameters of the superheat particles' 
initial distribution. So, let us specify the initial distribution in form: 

Af{r)^Ae-"^, (15) 

so that: 

oo 

= E ^ / ^nr)v^dv = ^Af (16) 
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is an initial conformal density of nonequlibrium particles' number. Let us cal- 
culate an initial conformal energy density of nonequilibrium particles: 



£i 



OO 

/ Aforr'dV = ^M. (17) 



Thus, we obtain for the initial average conformal energy: 

^0 = ^ = -. (18) 
AN a ^ ^ 

Let us calculate the parameter cto — eo/(£o+£i): where £o is an energy density of 
plasma's equilibrium component and £i is a nonequlibrium component's energy 
density. Plasma's equilibrium component is described via the distribution: 

fo = Bcxp = Bcxp I 1 = Bcxp ( — 1 . (19) 

Calculating conformal particles number densities and theirs energies relative to 
this distribution, we obtain: 

AAo = ^-^M^l^; io = (20) 



Comparing (j20p with an energy density of equilibrium plasma, which is deter- 
mined via its temperature by means of the relation (see [1]): 

£o=AA^, (21) 



we find: 



B = ^A/-. (22) 
Then: 

Using (HZl), dini) and ^ we obtain: 

A = {l-.o)^. (24) 

Then according to the referred above nonequilibrium particles' condition of 
smallness: 

ni{t) < no(t), (25) 
it has to be: 

a(l - ao) « ylN^^\ (26) 

On Figl7][9] the graphs of the initial staircase and exponential distributions at 
equal energy densities and particles' average energies for various nonequlibrium 
parameters are shown in comparison. 
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Figure 1: Comparison of staircase and exponential initial distributions at equal 
in both cases energy densities £i and average energies of particles Vq; a = 0.9. 
Staircase distribution is a thin line, exponential distribution is a heavy line. 
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Figure 2: Comparison of staircase and exponential initial distributions at equal 
in both cases energy densities ii and average energies of particles Vo; a = 0.5. 
Staircase distribution is a thin line, exponential distribution is a heavy line. 
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Figure 3: Comparison of staircase and exponential initial distributions at equal 
in both cases energy densities ii and average energies of particles Vq; a = 0.1. 
Staircase distribution is a thin line, exponential distribution is a heavy line. 



3 Research of nonequilibrium distribution's re- 
laxation 

Calculating function $(Z) relative to the distribution (IT5|) according to the 
formula (1131). we obtain: 



_ 9Z3 / 2{3Z + 6)Ko{2V3Z) 
iVSiS + 6Z)Ki{2V3Z)\ 



(27) 



27Z5/2 / ' 

where K^{z) is a Bessel function of the second sort (McDonald dunction) [S]: 



n{z) > 0,^{i^) > (28) 
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Relation of r and Z is determined via the following formula H 



(29) 



in which an obtained value of $(^) ([27| 

is to be substituted. On Fig. |3][7] the graphs showing an influence of initial 
distribution's parameters on plasma's temperature relaxation process are rep- 
resented. 
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Figure 4: Relaxation of plasma's temperature to the equilibrium: 
y = T{t)/Tq{t) subject to the parameter (Tq- bottom-up 
(To = 0, 01; 0, 1; 0, 3; 0, 5; 0, 9. Values of dimensionless time variable r are put 

on the abscissa axis. 
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Figure 5: Relaxation of plasma's temperature to the equilibrium: 
y = T(t)/Tq(t) for the parameter ctq = 0.9 for staircase function (dotted line) 
and Ae~"P function (firm line). Values of dimensionless time variable t are 



put on the abscissa axis. 
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Figure 6: Relaxation of plasma's temperature to the equilibrium: 
y = T{t)/Tq{t) for the parameter ctq = 0.3 for staircase function (dotted line) 
and Ae~"P function (firm line). Values of dimensionless time variable t are 

put on the abscissa axis. 
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Figure 7: Relaxation of plasma's temperature to the equilibrium: 
y = T{t) /To{t) for the parameter ctq = 0.01 for staircase function (dotted line) 
and Ae""*' function (firm line). Values of dimensionless time variable r are 

put on the abscissa axis. 

On Fig. |4] the results ol numerical modelling of equilibrium component's 
heating process by superheat particles for the initial distribution (jl5p in terms 
of equation (l29l) subject to the initial distribution's parameters are shown. It is 
obvious from this picture that temperature of plasma's equilibrium component is 
saturated up to the value To{t) at t ^ 3 — 4. On Fig I5]l7] temperature's relaxation 
process of the staircase and exponential distributions for various values of initial 
distribution's nonequilibrium parameter ctq is shown. It is clear from these 
pictures, that in the case of the exponential distribution temperature relaxes to 
the equilibrium rather slower than in the case of the staircase function, however 
qualitative behavior of 2/(r) functions's graphs coincides with each other. 

Using the results of equation (|29|) numerical integration in formula for the 
relaxation of the nonequilibrium distribution ([3|) , we obtain time dependence of 
distribution function's deviation from the equilibrium. On Fig|5] the results of 
equations' (|29l) and ([3]) simultaneous numerical integration subject to the value 
of initial distribution's nonequilibrium parameter are shown. 

Thus, an evolution of distribution in time variable Z does not depend explic- 
itly from the nonequilibrium parameter ao, but the relation of real time r and 
time variable Z do. This relation is given by formula (1291) and is represented in 
graphical form on Fig|8l 
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Figure 8: Relaxation of superheat component for distribution p5l) (A = 1) 
subject to time parameter Z: top-down - Z = 0.5, Z = 1, Z = 2, Z = A.Value 
of dimensionless time variable x is put on the abscissa axis, value lg(l + /) is 

put on the ordinate axis. 
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Figure 9: Dependence of distribution's maximum on dimensionless time t{Z), 
calculated by formula j^) relative to the initial distribution ijTSj) subject to 
initial distribution's nonequilibrium parameter cto- Top-down: ctq = 0.01, 
CTo = 0-1, ctq = 0-3, To = 5, ao = 0.9, 
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4 Conclusion 



Thus, carried out research has shown from the one hand though certain differ- 
ences the stabihty of nonequihbrium distribiition's relaxation scenario relative 
to parameters of the superheat particles' initial distribution and from the other 
hand has allowed to revcial the dc!peiidc;ncc; of initial distribution's average pa- 
rameters on Universe's evolution tinic!. As evident from aforecited pictures, 
maximum of nonequihbrium particles' energy spectrum shifts with time by the 
approximate law: 

T^max — • 

In conclusion. Authors express their gratitude to professors V.N.Melnikov, 
Yu.S.Vladimirov and D.V.Galtsov for the stimulating discussion of paper's mat- 
ters. 
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